Model for magnetic flux pattern induced by in-plane magnetic fields on spatially 
inhomogeneous superconducting interfaces of LaA103-SrTi03 
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The effect of spatial inhomogeneity on the properties of a two-dimensional non-centrosymmetric 
superconductor in an in-plane magnetic field is investigated, as it can be realized in LaA103-SrTi03 
interfaces. We demonstrate that the spatial variation of Rashba spin-orbit coupling (RSOC) yields 
a local magnetic flux pattern due to the field-induced inhomogeneous helical phase. For sufficiently 
strong fields, vortices can nucleate at inhomogeneities of the RSOC. 



The magneto-electric effects belong to the most in- 
triguing features of superconductors without inversion 
symmetry [U-0]- The connection between spin polariza- 
tion and supercurrents is one basic aspect which has nu- 
merous consequences. For example, spin polarization, 
by external magnetic fields or intrinsic ferromagnetic or- 
der, can induce a so-called helical phase in the supercon- 
ducting state. Such a helical phase can appear as sim- 
ple phase factor e lq r to the superconducting order pa- 
rameter, resembling a Fulde-Ferrell-Larkin-Ovchinnikov 
(FFLO) UlEl state 0-0, H, E3- However, also am- 
plitude modulations of the order parameters are possi- 
ble on this basis jfj EH Such magneto-electric ef- 
fects originate from spin-orbit coupling specific to non- 
centrosymmetric metals. 

Non-centrosymmetric superconductors (NCSC) most 
extensively investigated are those lacking mirror sym- 
metry with respect to a crystalline plane (e.g.: a-b- 
plane of a tetragonal crystal) and possess spin-orbit cou- 
pling of the Rashba-type with the well-known structure: 
ct(k x z) ■ S — agk • S Non-centrosymmetricity of this 
kind is realized in some heavy Fermion superconductors 
such as CePt 3 Si QJ] and CeTSi 3 (T=Rh, Ir) [U E3 and 
in superconducting thin films grown on a substrate. A 
particularly interesting example of the latter case is the 
superconductor found at interfaces between the band in- 
sulators LaAlC>3 and SrTiC>3 [l9j . Here, superconductiv- 
ity and spin-orbit coupling can be influenced by perpen- 
dicular electrical gate fields [2(J HH • This system provides 
a unique platform to study magncto-clcctric effects. Re- 
cently, scanning SQUID experiments have detected local 
magnetic patterns which have been interpreted as fer- 
romagnetic patches coexisting with the superconducting 
phase (HJ. Motivated by these experiments and the de- 
sire to detect signatures of the so far elusive helical SC 
phase we investigate magnetic properties of a 2D su- 
perconductor in an inhomogeneous environment exposed 
to an in-plane magnetic field or magnetization. 

A field parallel to the plane of a sufficiently thin film 
of a NCSC would generate an ideal helical phase without 
being disturbed by vortex lattice modulations. However, 
this phase would not leave obvious traces as its phase 
gradient can be compensated by a gauge field equivalent 
to the screening of magnetic fields in superconductors Q . 



The situation changes, if the NCSC is inhomogeneous in 
this plane, for instance, through modulation of spin-orbit 
coupling. We will show below on a simple geometry of 
two half-planes of different spin-orbit coupling strength, 
how in the helical phase magnetic flux pattern and even 
vortices oriented perpendicularly to the film can appear 
in the region of varying spin-orbit coupling (see Fig. [T] 
(d)). 

We model our system by the BCS-like Hamiltonian, 

n = *w(k)c ks ck, S ' - aJ2 st (q) s (q). 

k,s,s' q 

K ss , (k) = e k 5 ss < + <t ss > ■ [ag k + .g^sH] , (1) 

where cu_ s annihilates a quasiparticle state with momen- 
tum k and spin s, and 
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= 2 2^ H CT 2) ss 'C_ k+ s s c k+ s s >. (2) 

k,s,s' 

The single-electron spectrum is described by K ss i in- 
cluding the RSOC with strength a and the unit vector 
gk = kxz (k = k/|k|) and the Zeeman field j/jbH cou- 
pling to the spin {g: ^-factor and \xb'- Bohr magneton). 
The kinetic energy e k is measured from the Fermi energy 
e p . Finally, we restrict ourselves to a pairing interaction 
of strength a (> 0) in the (onsite) s-wave channel [24j|. 

The free energy derived from the above model can 
be expanded in powers of the gap function A (A(q) = 
a(_B(q))) to yield a generalized Ginzburg-Landau (GL) 
theory valid near T c [i-H i|. 

Tgl = J d 2 r [^ 2 |A| 2 +A 4 |A| 4 + iv|nA| 2 (3) 

(V x A) 2 



-K(z x H) • {ATIA + A(IIA)*} ■ 



where the coefficients are given by 

A 2 = N (in L + 2 7 .gV B # 2 



8tt 
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A 4 = A 7, K = N jv 2 F /2, and K = 5N gfi B v F j/2 with 
7 = 7C(3)/16(7rfc B T c ) 2 , N = (N+ + N-)/2 as the mean 
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density of states per spin, and 5N = N + — N- = 2aN^ 
(N±; density of states of the two bands split by the 

RSOC: = £k ± |agk| and Ng a measure for the 

particle-hole asymmetry of the Fermi surface). The co- 
variant gradient is given as II = — ihV + (2e/c)A, where 
A is the vector potential with VxA = 'Bint and Bi„j 
as the internal magnetic field in contrast to the Zee- 
man field H. The second term in A 2 describes the 
paramagnetic limiting effect through the Zeeman field 
H. The second gradient term involves magneto-electric 
effects arising from in-plane spin polarization. We ne- 
glect any spatial dependence perpendicular to the plane 
assuming the thickness of the SC interface in LaAlC>3- 
SrTiC>3 sufficiently small [l9l. [23j. Moreover, we restrict 
to \a\/Ep <C 1, and spatial variations of a(r) on scales 
larger than the Fermi wave length are incorporated by a 
straightforward generalization into the GL free energy. 

Before discussing the effects of these terms on the mag- 
netic properties of an inhomogeneous superconductor, we 
would like to briefly mention one important aspect of 
the magneto-electric effect and the corresponding helical 
phase. In a thin film with an in-plane magnetic field H 
paramagnetic limiting is the relevant pair-breaking mech- 
anism, as orbital pair-breaking is suppressed. The renor- 
malized upper critical field is given by 



H P (T) 



H P ,o(T) 



(5) 



y/1 -6N* /(AN*) 

with a gap function A(r) = Aoexp(iq • r), the he- 
lical vector defined as q = — (K/Kh)(z x H) — 
(2e/7ic)A, and the bare paramagnetic limiting field 

H p ,o = y/\ ln(T/T c )\/(giiB\/2rf). Obviously, the wave 
vector q is gauge dependent, whereas resulting measur- 
able quantities, such as the current densities j s + jz, de- 
fined below, are gauge invariant. The magneto-electric 
coupling stabilizes the superconducting state and in- 
creases the paramagnetic limiting field [f| Q. Conse- 
quently, the inhomogeneity of RSOC would yield inho- 
mogeneous nucleation of superconductivity in a uniform 
field even in systems which would show homogeneous su- 
perconductivity in zero field. 

A first idea of the magneto-electric effect can be gained 
by considering the GL equation obtained by variation of 
Tql with respect to A 



An .. 

V x B int = — (j s 
c 



" j-L ~t~ jme) 



(6) 



with current densities = i2eKTi{A* VA — c.c.} and 
]l = — (8e 2 K/c)\ A| 2 A, and the magneto-electric con- 
tribution j me = — 2eA|A| 2 (z x H). While j s and j L 
are standard currents due to phase gradient and Lon- 
don screening contributions, respectively, the last one is 
characteristic for magneto-electric phenomena due to the 
field- induced shifts of the Fermi surfaces jf| 0| . Note that 
for the uniform helical phase, the sum j s + j me + jz, van- 
ishes. The curl of Eq.© leads to the London equation 
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FIG. 1: (Color online) Spatial dependence of the Fourier com- 
ponents of the SC gap function c n (x) (a) and the vector po- 
tential a y , m (x) (b) along the x-axis in an inhomogeneous state 
[solid (red) curves] with and a homogeneous state [dotted 
(green) ones] without spatial modulation in the y direction at 
T/T c = 0.9, H/H Pt - = 0.976. (c) The free-energy densities of 
the two states are averaged over the y-range. The length are 
given in units of £(t) with d/£(t) = 3.0 and A (T)/£(T) = 9.0. 
(d) Geometry of the inhomogeneous 2D NCSC with strong 
(weak) RSOC on the right (left) half-plane whereby an in- 
plane magnetic field H|[x yields out-of- plane magnetic flux 
Bi„t|[l around the boundary (see Fig. [3]). 
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and A -2 = 32-7rA|A| 2 e 2 /c 2 as the London penetration 
depth. We find an effective source field whose first term 
originates from the ordinary phase gradient, while the 
second part is anomalous and depends on the spatial 
dependence of the Zeeman field, the spin-orbit coupling 
and/or the order parameter size [l6j . 

We now turn to the example of inhomogeneous spin- 
orbit coupling with two half-planes of different magnitude 
of 5 N described by 



SN (r) =SN 1 + r)( tanh(x/rf) + l) 



(9) 



with a boundary of thickness d. This gives rise to a differ- 
ence between K on the two sides: AK = g/j,BVF"fV^^o- 
The magnetic field H = (H, 0, 0) is uniform and di- 
rected along the x-axis (perpendicular to the domain 
boundary). We obtain the spatial structures of the 
SC gap function A and the vector potential A, by 
numerically solving the GL equations SFq^/SA = 
and (5J-"gl/^A = under suitable boundary conditions 
at x = +oo. Here, the Fourier transformations are 
defined by A(r) = fc_BT c J2 n c n(x) exp[z 2irny/Ly\ and 
A(r) = y(hc/2e)£(T)- 1 J2 m ay, m {x)exp[i2Trmy/L y ], 
and the boundary condition is expressed as 



i(+oo) = S„ 



kBT c 
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FIG. 2: (Color online) Spatial structures of the SC gap func- 
tion |A| 2 and the local current j = j s + jz, + j em (arrows) in 
the homogeneous helical state (a) and the inhomogeneous one 
with a spatial modulation in the y direction (b) (same param- 
eters as in Fig[TJ). The length are given in units of £(£). In 
(b), vortices threading the interface appear on the boundary 
at x = for the in-plane field H||x. 
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f(T) = hv FV / 'j/2\ \n(T/T c )\ as the coherence length, and 
ni, is defined as an integer minimizing |a a! o(+°o)|. Here, 
the parameters L y /£(T) = 115, SN /N = 0.16, and 
i] = 0.75 are used, and then, rib = —4 at T/T c = 0.9, 
H/H p ^ = 0.976 (L y : extension of the system along the 
y-direction) . 

We find two basic phases for the given geometry. One 
is a state homogeneous along y-direction with a single 
helical vector q y (+oo) and the other is an inhomogeneous 
state consisting of two main helical vectors q y (—oo) and 
q y (+oo). Figure[T]shows the numerically obtained spatial 
profiles of the Fourier components of A (a) and A y (b) 
along the x-axis in the homogeneous state [dotted (green) 
curves] and the inhomogeneous one [solid (red) ones] . 

In the homogeneous helical state, only one Fourier 
component c n whose wave vector corresponds to (7 y (+oo), 
ri = —4 = rib, is non- vanishing, so that |A| is uniform 
along the y-axis. The inhomogeneity of the RSOC is 
reflected not only in the x-dependence of |A| through in- 
homogeneous paramagnetic depairing effects (see Eq.([5])) 
but also in a large spatial variation of A y (x). The limit- 



ing difference in A y between the two sides is given by 

M^)-A v{ -oo) = ^H^ = ^H^ (11) 

where <&o is the standard flux quantum. The local cur- 
rent j = j s + ji + jem is shown in Figj2] (a). Near the 
domain boundary at x = 0, the current flows along the 
y-direction with a vanishing net current because of the 
Meissner screening effect. Consequently, as shown in Fig. 
[3] (a), the induced internal magnetic field is uniform along 
the y-axis, concentrated to the vicinity of the boundary 
on the length scale A and oriented along the z-axis. The 
out-of-plane magnetic flux located around the boundary 
is given by <fi — A y {-\-oo) — A y {— oo) per unit length of 
the boundary. 

The net magnetic flux is reduced for the other state 
when spatial dependence of |A| is introduced along the 
boundary (y-axis) . This means that several Fourier com- 
ponents in Eq. (|10[) are non-zero. For our choice of pa- 
rameters the gap function is a superposition mainly of 
the components c_2 and c_4 whose wave vectors n = — 2 
and n — —4 correspond to different helical vectors on 
the two sides q y (—oo) and ^(+00), respectively (FigfTJ). 
This reflects directly different RSOC and, thus, different 
shifts of the Fermi surfaces in the two half-planes @ . In 
the extreme case we expect here A y (+oo) — A y (— 00) = 0, 
such that 



Aq = q y (+oo) - q y (-oo) 



AK 



H. 



(12) 



Under this condition the two helical states have to be 
matched at the boundary which leads to a regular array 
of topological defects, corresponding to vortices lined up 
along the boundary. The spatial structures of | A | 2 and j 
are shown in Fig. [5](b). The vortices are visible through 
circular currents centered around the vortex cores where 
|A| is suppressed to zero. The internal field in this vor- 
tex state is shown in Fig. [3] (b). This vortex pattern is 
superposed on the magnetic field of opposite sign spread 
along the boundary. Note that in Eq.® now both con- 
tributions to B e /y |j z are active. 

The appearance of vortices corresponds here to a 
screening effect. The vortices being opposite to the back- 
ground magnetic field reduce the net magnetic flux on 
the boundary Naturally, at low in-plane fields ("weak" 
Fermi surface shift) the homogeneous helical state is sta- 
bilized and at a higher in-plane field a transition to the 
inhomogeneous state with vortices occurs. 

Insight on the relative stability of the two states can 
be gained through the free-energy density averaged along 
the y-axis shown in Fig. Q] (c) . Two prominent features 
can be identified, a rather localized contribution for the 
vortex phase (characterized through strong spatial varia- 
tions of the order parameter at the boundary) and a more 
extended contribution due to the magnetic flux which is 
more extended in the homogeneous phase. Energetically 
the transition between the two phases is determined by 
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FIG. 3: (Color online) Spatial distributions of the internal 
magnetic field given by Bj n t = z\/ x A y (x, y), in the homoge- 
neous state (a) and the inhomogeneous state with vortices (b), 
where Bi„t is normalized by (hc/2e)£ i (T)~ 2 and the length are 
given in units of 



the competition between the magnetic energy of the ho- 
mogeneous phase and the vortex energy. 

We provide here a simple estimate of the critical field 
for the transition between the two states. The energy 
of the homogeneous state is dominated by the mag- 
netic field contribution approximately given by Em 
(4>/2X) 2 (X/Att). The vortex state contributes mainly 
through the energy of each vortex e v . Per unit length 
along the boundary the energy corresponds to En) ~ 
e v /£ v , where £ v denotes the distance between the vortices 
and is derived through Aq — 2ir/£ v , i.e. each vortex in- 
cludes 2tt phase mismatch between the helical states of 
the two sides. Comparing these two energies yields the 
criterion 

H<S>lAK = 32ir 2 hKe v X. (13) 

With the independence of e„ and A through |A| 2 , we 
obtain the critical field for the discontinuous transition, 

tfc(T) = %g, f3 x = ^AKH p ^T), (14) 



where e v = \A\ 2 e v and A = X(\A(H = 0)|). Note that 
the spacing between the vortices is roughly proportional 
to the in-plane field H. 

The width d of the boundary has not appeared in our 
estimate, as we have assumed d <C X in our approxi- 
mation of E(ij . If d becomes comparable or even larger 
than A, the magnetic field distribution is spread out on 
the length d and E^\ ~ (cj) / 2d) 2 (d / 4n) such that for H c , 
X has basically to be replace by d in Eq. dTBl . Thus, the 
critical field increases. 

We address now some experimental aspects to detect 
the features discussed above. A direct way to observe the 
magnetic pattern introduced by the in-plane field is the 
local detection of out-of-plane magnetic fields (e.g. scan- 
ning SQUID microscope). Note that the sign of the mag- 
netic flux depends on the gradient of RSOC modulation. 
An island of stronger RSOC has boundaries of opposite 
magnetic flux patterns. No magnetic flux is generated 
in the region where the gradient of the RSOC is perpen- 
dicular to the applied in-plane field. A further signature 
of modulated RSOC is the inhomogeneous nucleation of 
superconductivity for an in-plane magnetic field which 
could be observed through resistivity measurements as a 
broadening of the superconducting transition. Moreover, 
we would like to address the critical field for the vortex 
state generation. Looking at the iJ-dependence of the 
free energy, we find that the in-plane (spin) magnetiza- 
tion would grow roughly linearly with H in the homoge- 
neous case and would discontinuously drop to a basically 
constant value for the vortex phase. 

Finally, we would like to comment that as one can see 
in Eq.®, not only inhomogeneous RSOC would gener- 
ate the flux pattern we discussed, but also other inhomo- 
geneities such as non-uniform in-plane fields would have 
a similar effect. Moreover, we would like to mention that 
twin-boundaries in non-centrosymmetric crystals could 
yields similar properties in a magnetic field [26]. In view 
of the fact that spontaneous ferromagnetism had been 
reported for the LaA103-SrTi03 interfaces, an in-plane 
magnetic field may not be necessary to generate flux pat- 
tern [12] • Generally, the configuration discussed here rep- 
resents a good way to detect features of the helical phase. 
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